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ABSTRACT 

The high energy continuum in Seyfert galaxies and galactic black hole candidates 
is likely to be produced by a thermal plasma. There are difficulties in understanding 
what can keep the plasma thermal, especially during fast variations of the emitted 
flux. Particle-particle collisions are too inefficient in hot and rarefied plasmas, and a 
faster process is called for. We show that cyclo-synchrotron absorption can be such a 
process: mildly relativistic electrons thcrmalize in a few synchrotron cooling times by 
emitting and absorbing cyclo-synchrotron photons. The resulting equilibrium func- 
tion is a Maxwellian at low energies, with a high energy tail when Compton cooling is 
important. Assuming that electrons emit completely self absorbed synchrotron radia- 
tion and at the same time Compton scatter their own cyclo-synchrotron radiation and 
ambient UV photons, we calculate the time dependent behavior of the electron distri- 
bution function, and the final radiation spectra. In some cases, the 2-10 keV spectra 
are found to be dominated by thermal synchrotron self-Compton process rather than 
by thermal Comptonization of UV disk radiation. 
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1 INTRODUCTION 

Particle-particle collisions among electrons, and electrons 
and ions, are rare in rarefied hot plasmas, such as the ones 
thought to be responsible for the high energy emission in 
Seyfert galaxies. It is therefore unclear if the plasma in these 
sources can be described by a thermal, Maxwellian distribu- 
tion, especially during the fast variations often observed in 
their X-ray fiux and spectrum. On the other hand the OSSE 
observations of the brightest Seyfert galaxies (i.e. NGC 4151 
and IC 4329, see e.g. Zdziarski, Johnson & Magdziarz 1996, 
Madejski et al. 1995) show the presence of a high energy 
cut-off, highly indicative of a thermal nature of the under- 
lying plasma. Similar indications are obtained from the sum 
of 60 combined spectra of 27 Seyfert galaxies (Zdziarski et 
al. 1995). 

In this paper, we propose a new thermalization mech- 
anism, based on the process of synchrotron and cyclotron 
self-absorption. Electrons embedded in a magnetic field emit 
and absorb their own cyclo-synchrotron radiation and can 
exchange energy by exchanging photons. We will show that 
complete thermalization occurs (in a few cooling times) in 
sources magnetically dominated and when the radiation is 



completely self-absorbed. Increasing the amount of Comp- 
ton coofing makes the peak of the Maxwellian distribution 
to shift to lower energies, with the development of a high 
energy tail. 

Using the found equilibrium electron distribution we 
compute the Comptonization spectra in the framework of 
the two-phase model of Haardt & Maraschi (1991) and 
Haardt, Maraschi & Ghisellini (1994). We calculate how the 
hot magnetized plasma, above a relatively cold accretion 
disk, reaches a steady state equilibrium distribution through 
synchrotron self-absorption and by Comptonizing the seed 
photons coming from the disk. We show that as long as 
the hot phase is magnetically dominated, the equilibrium 
temperatures and mean energies of the emitting plasma are 
in agreement with the existing observations, quite indepen- 
dently of the size of the region. There is however an impor- 
tant difference with respect to a pure thermal scenario, in 
which particles are assumed to always be in thermal equi- 
fibrium. Here, in fact we consider the case of continuous 
injection of new, energetic electrons, up to some maximum 
energy 7max»TicC^. This implies that: 

1) To reach steady state requires some form of escape 
of the injected electrons in order to avoid a pile up of cold 
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electrons. Also reacceleration could occur, but we do not 
consider this case here. 

2) The final distribution is sensitive to the mean energy 
of the injected electrons, i.e. to 7max, and to the spectral 
shape of the injected distribution. 

3) We neglect the possible role played by electron- 
positron pairs created by the produced 7-rays. This could be 
a self-regulating mechanism to keep 7max small (see Ghis- 
ellini, Haardt & Fabian, 1993). Here we assume a given fixed 

value of 7max. 

Note that the model proposed here is an hybrid be- 
tween non-thermal and thermal models, because it assumes 
a continuous injection of electrons (and not a re-heating), 
which nevertheless reach a thermal (or quasi-thermal) dis- 
tribution. 

Previous related studies on the subject of synchrotron 
reabsorption have been made by Ghisellini, Guilbert & 
Svensson (1988, hereafter GGS88) by Ghisellini & Svensson 
(1989), and by de Kool, Begelman & Sikora (1989). These 
studies demonstrated that a quasi-Maxwellian distribution 
can form at the low energy end of an otherwise power law 
distribution. In these cases, in fact, electrons were injected at 
very high energies (7max > 10^) and thin synchrotron emis- 
sion was the main cooling process. The problem of deriving 
the equilibrium distribution in the presence of self-Compton 
losses was briefiy discussed by Ghisellini & Svensson (1989). 

The present paper differs from previous work because 
we study a regime where cyclo-synchrotron emission is al- 
most completely self-absorbed, and the main cooling mech- 
anism is due to the Inverse Compton of photons produced 
in a relatively cold accretion disk. 

The structure of the paper is as follows: in Section 2 we 
describe the main assumptions of the model and the used 
equations, in Section 3 we describe our results, which are 
then discussed in Section 4. 



2 SET UP OF THE SYSTEM 

Assume that in a region of dimension R relativistic elec- 
trons are injected at a rate (5(7) [cm~^ s~^] between 7min 
and 7max. A tangled magnetic field B of energy density Ub 
makes them radiate synchrotron (S) photons. These pho- 
tons, together with photons produced externally to the re- 
gion, interact with the electrons by the Inverse Compton 
(IC) process. If the electron distribution extends to a 7max of 
the order of a few, the synchrotron spectrum is completely 
self absorbed, and the total radiation energy density {Ur) 
may be dominated by the photons produced externally. We 
will assume Ub 3> Ur- Due to self absorption and the pres- 
ence of external photons, this does not correspond to the 
prevalence of the synchrotron luminosity over the inverse 
Compton one. 

The electron distribution A''(7, t) is the result of the 
injection, the S and IC losses, and the energy gain due to 
self absorption (GGS88) 



7s and 7c are the cooling rates for synchrotron and Inverse 
Compton losses, respectively: 



7s 



4 (jTcp Ub. 



3 rrioC^ 



(2) 



(3) 



Equation (3) assumes that the Compton cooling is in the 
Thomson regime. The kinetic equation (1) difi'ers from the 
analougous equation (2) of GGS88 by the inclusion of the 
Compton loss term 7c. The factor H is defined below. The 
last term represents electron escape from the plasma region 
with the escape time being fcsc = R/Vesc- We use iicsc = c. 

We assume that the radiation energy density Ui is dom- 
inated by an external soft photon distribution, and by the 
'hard' radiation produced by Comptonizing these photons. 
The soft photons are assumed to arise from the reprocessing 
of half of the hard radiation by cold matter in the vicinity of 
the active region as detailed in Haardt & Maraschi (1991). In 
steady state, all the power injected in the form of relativis- 
tic electrons balances the escaping luminosity. Part of the 
power escapes as kinetic energy flux and part is converted 
into radiation, mainly at high (UV and X-ray) energies, via 
IC. The ratio of the two depends upon the ratio of the escape 
and the IC cooling timescales. 

Since the cyclo-synchrotron emission is completely self- 
absorbed, it is the IC process which is mainly responsible for 
the radiation losses, even if Ub S> Ur- The Compton process 
becomes efficient after a time « 2R/c, which is the timescale 
required to build up the soft photon radiation field. After 
this time the radiation energy density can be written as 



R 2 

Ur{t) = a— (l + TT)mcC 

c 



iV(7) 



Q(7)--^ (7-l)rf7, 

t-csc 

(4) 
where a is a numerical coefficient (of order unity) which 
depends on the geometry (~ 3/4 for a sphere and ~ 1/2 
for a slab). We use a = 3/4. The factor (1 + tt) accounts 
for the enhancement of the photon density due to Thomson 
scattering with the electrons, the optical depth of which is 



TT(i) ~ cttR 



Ar(7,t)d7 



(5a) 



At equilibrium, the number of injected electrons must be 
balanced by the number of escaping electrons. In steady 
state the optical depth becomes 



lim -rr(t) = aTRtasc 

t — »oo 



Q(7)d7 



(56) 



The factor H{-y,t) in equation (1) describes the heating 
of the electrons and their diffusion in energy. It is defined as 



H{j,t) = 



1 



J{u,t) 



3s{v,l)dv , 



(6) 
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where the momentum p is measured in units of rricC and 



where js (v, 7) is the power spectrum emitted by a single 
electron through the synchrotron process, and J{v, i) is the 
radiation mean intensity due to the S and IC process, with 
the former dominating the absorption (and hence the heat- 



© 0000 RAS, MNRAS 000, 000-000 



Thermalization in compact sources 



ing). J{p,t) is calculated by 



J{v,t) 






[1 



^-RK(,.,t)-l 



(7) 



where fi{i', t) and niy, t) are the emissivity and the absorp- 
tion coefficient of cyclo-synchrotron radiation, as calculated 
in GGS88. 

Since we are dealing with mildly relativistic electrons, 
we must specify the cyclo-synchrotron emissivity of the sin- 
gle electron. Most existing expressions are either nonrela- 
tivistic or only valid at large harmonics (Petrosian 1981). A 
simple one parameter phenomenological expression for the 
cyclo-synchrotron power spectrum is 



j^{v,-y) 



Ap utcUb. 



I'B 1 + 3p- 



■ exp 



2(1 - u/ub) 



l + 3p2 



(8) 



where vb is the Larmor frequency. This expression (i) has 
the correct cooling rate, equation (2), when integrated 
over frequency, (ii) has the correct frequency dependence, 
oc exp(— 2z^/t'_B), at large harmonics {f ^ ub) in the non- 
relativistic limit, (iii) has the correct frequency dependence, 
at large harmonics in the ultrarelativistic limit, and (iv) the 
agreement at smaller harmonics, i^/vb < 100, is better than 
40% at 7 = 2. 

For 7 > 2 we use the usual synchrotron formula av- 
eraged over isotropically distributed pitch angles (see, e.g., 
GGS88). 

The kinetic equation (1) can be solved for N{-y,t) as 
a function of (5(7), the magnetic field B, and the size R. 
Equivalently, besides R, one can specify the injected and 
the magnetic compactnesses. The total (including rest mass) 
injected compactness is 



where 



f- ■ — 

•tmj — 
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(10) 



and where V is the volume of the region. The magnetic 
compactness is defined as 



4 






RUb 



(11) 



In the two phase model for Seyfert galaxies, the high 
energy radiation could be produced by a uniform corona 
(Haardt & Maraschi 1991) or by localized active regions 
(possibly highly magnetized; Haardt, Maraschi & Ghisellini, 
1994). Assuming the second possibility, the dimensions of 
each region, or blob, is of the order of a few Schwarzschild 
radu. (i.e. 7? ~ 10^^ cm). 

Observations of the high energy continuum and vari- 
ability in Seyfert galaxies indicate that their compactness 
is of the order of 1-300 (Done & Fabian, 1989). We then 
assume that ^inj is in these range of values. Our assumption 
of Ub 3> Ur then also fixes the range of Ib- When escape is 
not important (i.e., for i?B ^ 1, see below) we have 



Ub 



167r Ib 



1- 



1 



< 7 > 
where < 7 > is the average of 7 over Q("f) 



(1 + tt) 



(12) 



The ratio of the cooling time scale, fcooi 
7c), to the escape time, fesc, is given by 

4~r 



t'cool 

t'CSC 



(l+7)^B(l + f/./l7B) 



(7-l)/(7s + 



(13) 



This is always less than unity for £b ~ 1, implying that elec- 
trons cool rather than escape for sufficiently large compact- 
nesses. Below we use Ib = 10 or 30. 

Electrons may also thermalize through Coulomb en- 
ergy exchange (e.g.. Stepney 1983, Dermer & Liang 1989, 
Nayakshin & Melia 1996, Pilla & Shaham 1997). The ther- 
mal electron-electron Coulomb energy exchange rate can be 
approximated with (e.g.. Stepney 1983) 



7coui(6) 



ttcIu A 



8i?ei/2(vri/2 + ei/2) 



(14) 



where In A is the Coulomb logarithm, typically ~ 10-20, and 
O = kT/nieC? is the dimensionless temperature. The ther- 
mal average of equation (2) becomes 

7s(e) ^ (c/i?)(4/3)e(i + 4e)fe . (i5) 

Thermalization by synchrotron selfabsorption dominates 
when 7s > 7coui, which gives Ib > (9/327r^''^)TT In A/O^''^ 
for O < 1. One sees that the Coulomb process dominates for 
small temperatures and large tt (i.e., large electron densi- 
ties). From equations (5b), (9), and (10), we find that the 
Thomson optical depth is given by 



TT = 



3 Vc: 



An 



< 7 > 



(16) 



Thermalization by synchrotron selfabsorption then domi- 
nates for temperatures 



e >0.11 



In A lin 

< 7 > -^B 



2/3 



(17) 



3 RESULTS 



3.1 Time evolution 



Figure 1 shows the evolution of the electron distribution to- 
wards the equilibrium configuration. The injected electrons 
have a Gaussian energy distribution peaking at 7 = 10. The 
magnetic field has a value B ~ 5.5 x 10'^ G and the in- 
jected power corresponds to £inj = 1. The size of the source 
is i? = 10^^ cm. As can be seen, the equilibrium shape of 
r(p) = aTRN{p) is reached in a time ~ O.lR/c, about equal 
to the nonrelativistic cooling time (see eq. 13). With the as- 
sumed input parameters, the synchrotron terms (emission, 
absorption and energy diffusion) in the kinetic equation are 
dominant over Compton losses. Gains and losses in this case 
almost perfectly balance. As a result the equilibrium electron 
distribution is a Maxwellian. Figure 1 also shows that the 
high energy part of the Maxwellian distribution is formed 
earlier than the low energy part, due to the higher efficiency 
of exchanging photons of the high energy electrons. A slower 
evolution takes place after 0.1(R/c), as balance between elec- 
tron injection and electron escape is achieved on a time scale 
of a few fosc. Only then have both the shape and the ampli- 
tude of N{p) reached their equilibrium values. 
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4=30, /i„, = 0.1, 0.3, 1, 3, 10, 30, 100 
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Figure 1. Evolving electron distribution at different times (mea- 
sured in units of R/c), as labelled. The size of the source is 
R = 10^^ cm, ^B = 10, and £inj = 1. The injected distribution is 
a Gaussian centered at 7 = 10. 



Figure 2. Equilibrium electron distributions for different injected 
compactnesses (decreasing from top to bottom). R = 10^^ cm and 
€b = 30 are assumed. 



3.2 Steady equilibrium distributions 

The equilibrium distributions for different values of the in- 
jected power are shown in Figure 2. The magnetic compact- 
ness is set to (.B. = 30, corresponding to B = 9.6 x 10^ 
G for R = 10^^ cm. In all cases, the injected distribution 
is (3(7) = Qop/7 exp(— 7/7c), where Qo is a normaliza- 
tion constant and 7c = 3.33. The minimum Lorentz factor 
is 7min = 1- The resulting mean injected Lorentz factor is 
< 7 >~ 4.6. As long as ^inj <C ^b, the distribution is a 
quasi-perfect Maxwellian at all energies. This is the result 
of the quasi-perfect balance between energy gains and losses, 
as also seen in Fig. 2: Compton losses (not balanced by a 
corresponding Compton heating term) are a small perturba- 
tion. As ^inj increases towards values ~ ^b, Compton losses 
start to be relevant, competing with synchrotron processes. 
At high energies, losses overcome gains, and the electrons 
diffuse backwards in energy, until the subrelativistic regime 
is reached. Only in this regime, the increased efficiency of 
synchrotron gains (relative to losses) halts the systematic 
backward diffusion in energy, and a Maxwellian can form 
(see Ghisellini & Svensson 1989). We can evaluate the tem- 
perature of this part of A'^(7) [or, equivalently, of N[p)] by 
fitting a Maxwellian to the the low energy part of the dis- 
tribution, up to energies just above the peak. 

The resulting temperatures are plotted in Figure 3 as 
a function of ^inj. As can be seen, for low values of £inj the 
temperature is quasi-constant, and decreases for li-nj greater 
than 1. Also shown in Figure 3 is the behavior of the "effec- 
tive temperature" Toff, defined by [Ocft = kTcS / {rnc(?)\. 

4 



leGeff -f 460 



< 7 - 1 >, 



(18) 



where the average is over the electron distibution N{j). As 
expected, for small ^inj, TcS is well approximated by T, while 
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Figure 3. Temperatures and effective temperatures of the dis- 
tributions shown in Fig. 2. Also shown are the corresponding 
values of the Comptonization y parameter. Above the solid line, 
synchrotron self absorption dominates over Coulomb exchange as 
the thermalization mechanism. 



the effective temperature is larger than T for ^inj ^ 10, as a 
result of the shift of the Maxwellian part to small energies, 
and the corresponding appearence of the high energy tail. 
Along the sequence, the optical depth increases from 

TT = 5 X 10"^ for ^inj = 0.1 to TT = 5 for ^inj = 100, 

according to equation (16). For our parameters and InA = 
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20, equation (17) becomes Q > 0.03(^inj)^''^, which is plotted 
as the solid line in Figure 3. One sees that that synchrotron 
self absorption dominates the thermalization for all cases 
with linj smaller than about 60 . For the case linj ~ 100, one 
cannot neglect Coulomb thermalization. 

With the values of Oeff and rx, it is of interest to cal- 
culate the Compton parameter y defined by 



y = (rT + r^)(16e,\+4eeff) 



(19) 



which is also plotted in Figure 3. At equilibrium, tt is given 
by equation (16) and 0cff cannot exceed the value of 0cft 
obtained by doing the average in equation (18) over the in- 
jected distribution Q{'j). This implies that y has a max- 
imum possible value, j/max- This limiting value is reached 
only when radiative cooling is negligible. The actual value 
of y is then further constrained by the amount of radiative 
cooUng. 

For .^inj > 1, j/max S> 1, but the amount of Compton 
cooling (dominant in this cases) keeps the value of y close 
to unity. The fact that y ~ 1 is to be expected: this value 
ensures equal power between the Comptonized radiation and 
the soft seed photon emission. 

For £inj < 1, ymax is instead less than unity: Compton 
coohng is negligible, and the actual value of y is only slightly 
less than j/max, due to synchrotron losses, which, albeit small, 
are in these cases dominant over Compton losses. 



3.3 Spectra 

In Figure 4, we show the radiation spectra corresponding to 
some of the equilibrium electron distributions of Figure 2. 
Each spectrum consists of several continuum components: 

i) the self-absorbed synchrotron spectrum (S); 

ii) the Comptonized synchrotron spectrum (SSC); 

iii) the thermal soft component (bump); 

iv) the spectrum resulting from Comptonization of the 
photons from the thermal bump (IC); 

v) the Compton reflection component. 

In addition, not shown in the plots, fluorescence line 
emission from photoionized iron in the accretion disk and 
absorption features due to partially ionized warm gas along 
the line of sight are expected to be present in the spectrum. 

For each model the values of the compactnesses of the 
Compton (SSC+IC) component and the synchrotron com- 
ponent are calculated as discussed in §2. In the adopted ge- 
ometry, neglecting the disk albedo, the compactness corre- 
sponding to the soft photon input is ^aoft ~ ^ava/(l + e"^'"'), 
where ^ava is the compactness available for radiation pro- 
duction, defined as the difference between the injected and 
the escaping kinetic compactnesses, i.e. subtracting the rest 
mass energy of the electrons from the total energy. 

The temperature of the soft radiation is then calculated 
assuming black body emission from a region of radius R. The 
synchrotron and soft component are then Comptonized in a 
plane parallel slab following the prescriptions described in 
Haardt (1994). For the synchrotron component we assume a 
homogeneous source function within the scattering medium, 
while the thermal soft photon input is localized below the 
active region, giving rise to anisotropic Compton emission 
(Haardt 1993). A face-on line of sight is assumed for all the 
spectra. 

The following features may be noticed: 



For linj ^ 1, the Comptonized spectrum is bumpy, due 
to the small value of tt and the large value of O. As re- 
ported in Figure 3, the y parameter is smaller than unity, 
making Compton losses relatively unimportant. The syn- 
chrotron component, albeit self-absorbed, is more impor- 
tant than the Comptonized power. The luminosity in the 
thermal "bump" (in the UV) is roughly half of the syn- 
chrotron luminosity (in the IR). Therefore, SSC and the IC 
components have approximately the same luminosity ratio 
as between the synchrotron and the thermal "bump" lumi- 
nosities. |j Due to the anisotropy Compton effect, the 2-10 
keV band is dominated by the SSC component, rather than 
by the IC. This is in contrast to the general view interpret- 
ing the high energy emission in Seyfert galaxies as due to 
Comptonization of thermal bump photons. Here the thermal 
bump and the X-ray flux are not directly related. 

Increasing ^inj up to 1, the Compton cooling becomes 
more important {y approaches 1), making i) the the syn- 
chrotron component (and the related SSC) weaker, and ii) 
the hard X-ray slope flatter. 

For linj ^ 3 the Comptonized spectra have the typical 
IC power law shape, with increasing (steeping) spectral in- 
dices as linj increases. The reason for this steepening even 
for constant j/ = 1, is explained in detail in Haardt, Maraschi 
& Ghisellini (1997). 

Finally, note that for linj ^ 3, the high energy part can 
be described by an exponential, since the electron distri- 
bution is a quasi-perfect Maxwellian in the entire energy 
range. For linj ~ 3, the electron distribution is more complex 
(see Fig. 2), resulting in a more complex radiation spectrum 
above ~500 keV. 



3.4 Stellar compact sources 

The relevant quantities for the synchrotron thermalization 
process are the injected and the magnetic compactnesses. 
These parameters are independent of radius, and therefore 
the process should be largely independent upon the size of 
the source. However, since B oc (Ib/R)^ (equation 11), 
the same value of Ib now implies B = 9.6 x 10^ Gauss. 
In turn, the increase of the magnetic held implies that the 
energy range of the electrons emitting optically thick syn- 
chrotron radiation is reduced. This affects the synchrotron 
emission for small values of linj , which now can have an opti- 
cally thin component from the most energetic electrons. For 
illustration, we have repeated our calculations assuming a 
source size of lO'^ cm, i.e. few Schwarzchild radii for a galac- 
tic black hole. For linj ~ 1 the resulting equilibrium electron 
distributions are almost identical to the R = 10^^ case, since 
the synchrotron radiation is completely self-absorbed. For 
linj ~ 1, the results are similar, but the resulting effective 



* Note that in Fig. 4 the thermal "bump" appears instead as lu- 
minous as the synchrotron component: this is due to the assumed 
isotropy of the synchrotron emission, while the "bump" emissiv- 
ity has a cos i angular pattern, where i is the viewing angle (here, 
1 = 0°). Note also that, while the SSC emission is isotropic at all 
energies, the IC component is not: the emission of the first or- 
der scattering in the i = 0° direction is depressed (Haardt 1993). 
At higher energies, both the IC and the SSC components are 
isotropic, and here their luminosity ratio has the expected value. 
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Figure 4. Radiation spectra calculated with the electron distributions shown in Fig. 2. Dashed line: reprocessed thermal bump and IC 
components. Dash— dotted line: synchrotron and SSC components. Dotted line: Compton reflection component. Solid line: total spectrum. 



temperatures are smaller (by a factor 2 for the li^j =0.1 
case), due to the presence of a thin synchrotron component. 
The computed radiation spectra are shown in Figure 
5. They exhibit the same basic features as seen in the ex- 
tragalactic case (Fig. 4), but with two major differences: i) 
as expected, the synchrotron emission is more important, 
with a corresponding more relevant SSC component con- 
tributing to the X-ray spectrum for l^nj ^ 3. The increased 
synchrotron cooling makes the Comptonized components 
steeper than the corresponding cases in Figure 4; ii) the 
peak energies of the thermal "bump" and the synchrotron 
emission are now in the soft X-ray and UV bands, respec- 
tively. 



4 DISCUSSION AND CONCLUSIONS 

We have shown that the exchange of cyclo-synchrotron pho- 
tons in a self absorbed source is a very efficient thermalizion 
mechanism, especially when the density and temperature of 



the plasma prevent the thermalization via particle-particle 
processes. 

In our calculations we assume that the magnetic energy 
density in the flaring corona dominates over both gas and 
radiation pressure. This is sensible assumptions in the pic- 
ture of an active corona. In fact the accretion mechanism 
is thought to amplify the magnetic field via some sort of 
dynamo process, and eventually particles in the corona are 
energized via buoyancy and reconnection of the field lines. 
Regardless of the largely unknown physics of such processes, 
it is conceivable that the particles, (and the radiation field 
they create because of cooling), can not be energetically 
dominant, as their energy is ultimately drained from the 
(amplifield) magnetic field. Note that this argument is true 
only if the reconnection timescale is larger than the cooling 
timescale for a radiating electron, something we have im- 
plicitly assumed allowing for steady state particle injection. 

We now discuss how the main results of the model 
change by varying some of the input parameters. 

As already mentioned, we have assumed, for all 
shown cases, an injected electron distribution Q{-f) oc 
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Figure 5. Radiation spectra for the same values of £b and £inj as in Fig. 4, but for i? = 10^ cm. Dashed line; reprocessed thermal bump 
and IC components. Dash-dotted line: synchrotron and SSC components. Dotted line: Compton reflection component. Solid line: total 
spectrum. 



xp/7 exp(— 7/7c) with 7c = 3.33, resulting in < 7 >~ 4.6. 
The emitted radiation spectrum is basically controlled by rr, 
since its value, coupled to the condition j/ ~ 1, determines 
the physical conditions in the source. Since tt oc ^inj/ < 7 > 
(eq. 16) , the assumed shape of the injected distribution is not 
important, as long as different distributions have the same 
< 7 >. Changing £inj and < 7 > by the same factor hence 
gives rise to the same photon spectrum, provided that the 
condition Is 3> ^inj is matched, and that the synchrotron 
emission is mostly self-absorbed. In other words, the same 
spectra shown in, e.g., Figure 4, would be obtained for values 
of £inj different than those indicated if an input distribution 
with a different < 7 > were used. It is important to note 
that the value of < 7 > cannot be taken too large, as it 
changes the fraction of self-absorbed synchrotron radiation. 
As long as the synchrotron emission is mostly self-absorbed, 
the spectra obtained for the same value of tt and ^inj/ < 7 > 
are indeed identical. Going to very large < 7 >, the condi- 
tions of self-absorption breaks down: the electrons cool down 
mainly because of synchrotron losses, with neglibible Comp- 



tonization, changing completely the resulting radiation spec- 
trum. 

As can be inferred from the shown cases, the value of 
R is largely unimportant: its value determines the peak fre- 
quencies of the synchrotron and bump components, but not 
the main characteristics of the electron and radiation spec- 
tra. 

The thermalization process described in this paper is 
operating whenever energetic electrons and some magnetic 
field is present, but it becomes of great interest in the study 
of the high energy spectrum of compact sources when the 
mean energy of the electrons is less than a few MeV and 
when the magnetic field is energetically dominant. These 
conditions ensure that the synchrotron spectrum is com- 
pletely self-absorbed, and the entire electron distribution is 
influenced by self-absorption. This also ensures that the syn- 
chrotron emission does not overproduce the observed IR (op- 
tical) emission observed in radio-quiet AGN (galactic black 
hole candidates). However, the cyclo-synchrotron process 
does contribute to the IR (optical) emission, thought to be 
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mainly due to reprocessing of the primary disk radiation by, 
e.g., dust. 

Since the same electrons contribute also to the X-ray 
band, simultaneous variability in the IR and X-ray bands 
should clearly identify this component, allowing the deter- 
mination of the fraction of synchrotron radiation. Its relative 
importance decreases by increasing the injected compactness 
(see Fig. 4). 

An important outcome of our calculations concerns the 
X-ray emission. If the magnetic field is dominant (i.e. its en- 
ergy density is larger than the radiation energy density) and 
the Comptonization parameter y is smaller than unity, the 
scattering of internal cyclo-synchrotron photons is more im- 
portant than the scattering of the external, "disk" photons. 
The variability pattern expected in our model is therefore 
rather complex, since different components can contribute 
in different bands, and their relative importance depends on 
the magnetic field dominance. When synchrotron emission 
is important, the 2-10 keV emission could be only indirectly 
related to the UV bump photons. Simultaneous UV and X- 
ray variability (though with different amplitudes) should al- 
low the determination of the importance of the synchrotron 
and SSC components. 
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